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摘 要 : 内 蒙古 乌拉 羡 河 流域 是 典型 的 干旱 半 干 旱 牧 


区 草原 内 陆 河流 域 ,生态 系统 极为 脆弱 ,气候 变 暖 和 人 类 活 


动能 直接 影响 水 文 变化 。 本 研究 采用 改进 的 SWAT 水 文 模型 .M-K 趋 势 检 验 .降水 -径流 量 双 累积 曲线 以 及 情景 分 


析 等 方法 ,系统 分 析 了 1981 一 2020 年 乌拉 盖 河 流域 径流 时 空 变 化 特征 ,并 量化 流域 不 同时 期 .不 同 河 段 气候 变化 和 


人 类 活动 对 径流 影响 的 差异 。 结 果 表 明 :SWAT 模 型 在 乌拉 盖 河 流域 的 适用 性 良好 , 率 定 期 及 验证 其 
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在 0.62 以 上 ,PBLAS 小 于 18.8%。 近 40 a 在 流域 呈 暖 干 化 趋势 下 ,径流 量 在 上 中、 下游 均 呈 显著 减少 趋势 ,上 且 在 
2000 年 发 生 突变 。 气 候 变 化 和 过 度 放牧 、 盲 目 开 坚 以 及 水 利水 库 建 设 等 人 类 活动 对 流域 径流 变化 的 贡献 率 分 别 为 
95.84% 和 4.16%。 人 类 活动 对 流域 不 同 河 段 的 贡献 率 也 有 差异 ,由 上 游 到 下 游 的 贡献 率 分 别 为 1.69% .4.36% 和 
5.03%。 且 在 不 同时 间 段 内 的 贡献 率 也 有 较 大 的 差异 ,由 1980 年 的 88.26% 减 少 到 2020 年 的 25.47% ,不 同时 间 段 内 
不 同人 类 活动 方式 导致 径流 变化 趋势 及 幅度 也 有 差异 。 人 研究 结果 能 为 牧区 草原 内 陆 河 流域 水 资源 的 可 持续 利用 


及 合理 调度 提供 参考 依据 。 


关键 词 : 径流 ; 气候 变化 人 类 活动 ; SWAT 模 型 ;乌拉 盖 河 流域 


近年 来 ,受气 候 变 化 及 人 类 活动 加 剧 的 双重 影 
啊 , 水 资源 在 时 空 上 重新 分 配 且 洪 小 ,干旱 等 灾害 
频 发 ,特别 是 在 干旱 、 闪 干旱 区 草原 内 陆 河 流域 万 
为 明显 一 。 水 资源 作为 制约 草原 内 陆 河 流域 发 展 
的 重要 资源 ”” ,其 变化 无 疑 对 当地 生态 系统 安全 及 
可 持续 发 展 产 生 重 大 影响 。 因 此 ,深入 剖析 流域 水 
文 循环 的 特点 ,厘定 气候 变化 和 人 类 活动 对 径流 的 
影响 ,能 够 为 最 大 限度 减少 环境 变化 对 径流 产生 的 
不 利 影响 提供 理论 依据 ,对 干旱 区 草原 内 陆 河 流域 
的 可 持续 发 展 具有 重要 的 现实 意义 "。 

目前 定量 区 分 气候 变化 和 人 类 活动 对 径流 贡 
献 的 方法 主要 有 统计 分 析 法 .水量 平衡 法 和 水 文 模 
型 模拟 法 ”。 统 计 分 析 法 依赖 于 大 量 的 历史 实测 数 
据 ,因此 ,在 缺乏 水 文 资料 的 地 区 难以 使 用 。 而 水 
量 平 衡 法 主要 根据 水 量 平 衡 建 立 数学 模型 ,通过 改 
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变 模型 中 的 单个 因子 来 估算 各 因子 的 水 文 效应 ; 
此 ,该 方法 无 法 了 解 水 文 循环 物理 过 程 "H。 水 文 模 
型 因 其 能 考虑 流域 空间 异 质 性 且 对 水 文 循环 物理 
过 程 进行 刻画 而 成 为 目前 最 佳 。 如 刘 柏 君 AY 
4°) Aawar 采用 水 文 模拟 途径 ,定量 评估 了 气候 
变化 和 人 类 活动 对 不 同 流域 径流 变化 的 影响 。 然 
而 ,目前 有 关 径 流 演变 特征 及 模拟 预测 的 研究 多 数 
集中 在 世界 较 大 的 河流 及 其 支流 ,很 少 有 学 者 关注 
草原 内 陆 河流 域 ,并 极 少 对 此 类 流域 从 定量 的 角度 
进行 研究 。 

乌拉 盖 河 流域 作为 内 蒙古 最 大 的 内 陆 河 流域 ， 
是 乌 珠 穆 沁 草原 的 重要 组 成 部 分 ,同时 也 是 内 蒙古 
自治 区 重要 的 畜牧 业 基 地 呈 ,不 仅 为 该 地 提供 可 贵 
的 水 资源 ,而 且 发 挥 较 大 的 生态 效益 ,是 乌 珠 穆 沁 
草原 的 生命 线 , 也 是 维系 草原 生产 力 的 重要 生态 因 
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素 ”"*“。 此 外 ,由 于 气候 变化 及 人 类 活动 的 持续 影 
啊 ,流域 草原 面积 不 断 下 降 ,草原 沙化 以 及 盐 碱 化 
的 程度 逐渐 加 重 , 对 牧区 草原 生态 平衡 及 经 济 发 展 
带 来 了 极 大 的 威胁 “。 然 而 ,目前 针对 乌拉 盖 河 
流域 径流 变化 及 其 主要 驱动 要 素 的 分 析 鲜 有 报道 ， 
特别 是 2000 年 实施 水 土 保 持 政策 以 来 流域 径流 变 
化 特征 及 其 原因 尚 不 清晰 。 因 此 ,本 研究 选取 乌拉 
盖 河 流域 为 研究 区 ,结合 改进 的 SWAT 水 文 模型 ""、 
情景 设置 ”及 数理 统计 等 方法 探讨 乌拉 盖 河 流域 
径流 的 时 空 变化 规律 ,定量 评估 不 同时 段 气 候 变 化 
和 人 类 活动 对 流域 径流 的 影响 。 人 研究 结果 能 够 深 
化 对 干旱 、 半 干旱 草原 内 陆 河 流域 生态 水 文 过程 及 
演变 规律 的 认识 ,为 内 陆 河 流域 水 资源 的 可 持续 利 
用 及 合理 调度 提供 参考 依据 。 


1 研究 区 概况 


乌拉 善 河流 域 (116?20'~11959 下 ,44*02'~46"42'N ) 
位 于 乌 珠 称 沈 草原 ,海拔 796~1937 m( 图 1)。 乌 拉 
盖 河 流域 面积 为 3x10' km ,是 内 蒙古 最 大 的 内 流 河 
流域 ,河道 长 度 为 2x10; km, 多 年 平均 径流 量 ; 
4.17 m°- s~ (1981—2012 年) ,并 且 在 河道 下 游 孕 
育 了 乌拉 盖 湿 地 2 。 湿 地 地 下 水 资源 是 人 畜 用 水 、 


小 型 分 散 的 农田 种 植 灌溉 用 水 ,城镇 用 水 及 部 分 工 
矿业 用 水 的 主要 来 源 ”。 乌 拉 盖 河流 域 属于 半 湿 
润 . 半 干旱 大 陆 性 气候 ,年 平均 降水 量 由 东 向 西 递 
减 , 约 为 250~400 mm, 年 平均 温度 约 为 -0.9 Co), 
流域 内 的 地 带 性 植被 主要 为 典型 草原 ,流域 东北 部 
为 草 甸 草原 植被 ,西南 部 为 典型 草原 ,中 部 为 两 者 
过 渡 地 带 ' 中 。 为 了 分 析 马 拉 羡 河流 域 不 同 河 段 的 
水 文 情势 ,将 河道 分 为 上 游 中游 和 下 游 ,乌拉 盖 水 
库 以 上 为 上 游 ,由 乌拉 羡 水 库 到 胡 稍 庙 为 中 游 、 胡 
稍 庙 到 索 林 淖尔 为 下 游 ”“。 


2 数据 来 源 与 研究 方法 


2.1 数据 来 源 

为 了 构建 乌拉 盖 河 流域 水 文 模型 ,采用 的 基础 
数据 包括 数字 高 程 模型 (DEM) .土地 利用 / 覆 被 £ 
BOSH .气象 数据 及 实测 径流 数据 。DEM 数据 来 源 
于 美国 国家 航空 航天 局 (http://www.nasa.gov) ,空间 
分 辨 率 为 30 m。 为 了 有 效 揭示 地 表 和 覆盖 变化 的 影 
响 ,如 过 度 放牧 造成 的 草地 退化 等 ,土地 利用 / 履 被 
数据 采用 二 级 分 类 ,分 别 为 林地 、 高 覆盖 度 草地 、 中 
覆盖 度 草地 、 低 覆盖 度 草 地 水 域 . 耕 地 、 建 设 用 地 


和 沙 地 8 类 ;空间 分 辩 率 为 1 km, 由 中 国 科学 院 环境 
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图 1 乌拉 盖 河 流域 水 系 及 站 点 分 布 


Fig. 1 Distribution map of water System and station in the Ulagai River Basin 
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资源 与 数据 中 心 免 费 下 载 (http:/www.resdc.cn) , 包 
括 1980 年 .1990 年 .2000 年 .2010 年 和 2020 年 共 5 
期 。 土 壤 数 据 从 世界 土壤 数据 库 (http://westdc.west- 
gis.ac.cn/data/) 获 取 , 空 间 分 辨 率 1 km。 气 象 数 据 采 
用 中 国 气象 数据 网 (http:/data.cma.cn) 提供 的 
1981 一 2020 年 乌拉 盖 以 及 周边 锡林浩特 市 .和 东 乌 珠 
穆 沁 旗 等 12 个 气象 站 点 的 逐日 降水 .最 高 /最 低 气 
温 、 风 速 .相对 湿度 等 数据 。 由 于 牧区 草原 为 缺 资 
料 区 域 ,站 点 实测 数据 稀少 。 因 此 验证 模型 时 采用 
实测 径流 数据 .实地 采样 数据 .MODIS 积 雪 数据 及 
蒸 散发 数据 进行 多 目标 率 定 。 实 测 径流 数据 由 
1981 一 2003 年 奴 乃 庙 水 文 站 和 2004 一 2012 年 乌拉 
盖 水 库 人 库 水 文 站 合成 的 逐 月 径流 量 资料 。 采 样 
数据 为 2019 一 2020 年 实测 水 位 .断面 等 换算 成 为 径 
流量 。 积 雪 数 据 和 蒸 散 发 数据 分 别 选择 2000 一 
2012 年 MODI 的 MODI10A1 和 MOD16A2 产 品 , 其 时 
空 分 辩 率 均 为 8 d,500 m。 通 过 美国 NASA 陆地 过 
程 分 布 式 数据 档案 (https:Wlpdaacsvc.cr.usgs.gov/ap- 
peears/) 免 费 下 载 。 

2.2 SWAT 模 型 及 模拟 精度 评定 

由 美国 农业 部 和 农业 研究 局 开发 的 半分 布 式 
(Soil and Water Assessment Tool ,SWAT) 水 文 模型 能 
很 好 的 模拟 流域 水 文 循环 过 程 并 能 进一步 量化 流 
域 水 文 过 程 对 变化 环境 的 响应 i。 本 人 研究 采用 
队 自 主 改进 的 SWAT 模 型 ,包括 改进 的 融雪 模块 和 
子 流域 分 割 法 ”, 通 过 增加 累积 温度 确定 条 件 来 
分 离 降雨 和 降雪 类 型 以 及 在 流域 划分 时 增加 土地 
利用 / 歼 被 变化 节点 来 提高 模型 模拟 精度 。 

为 了 提高 模型 的 适用 性 ,以 水 文 站 实测 数据 为 
基础 ,应 用 SWAT-CUP 软 件 的 SUFI-2 算 法 进行 模型 
参数 敏感 性 和 不 确定 性 分 析 ”” 。 选 取 纳什 效率 系 
数 (NES) , 均 方 根 误差 (PBLAS) FE ABR) GA 
个 评价 指标 对 模型 模拟 径流 结果 进行 评价 ”。 此 
外 ,以 MODIS 积 雪 面积 和 薰 散 发 数据 对 模型 模拟 的 
职 雪 面积 和 蒸 散发 量 曲线 的 拟 合 度 进 行 模拟 结 
的 多 源 验 证 ,以 保证 模型 的 适用 性 及 稳定 性 。 评 价 

旧 标 公式 如 下 : 


NSE=1- 


L” -QY 
yo," = prey 


zomg") 
x0 


(1) 


PBIAS= x 100 (2) 
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OO sorso] 


pze- hze- Ee] 
(3) 

式 中 : 0,” 为 实测 值 (ms ); 0 为 模拟 值 (mi…s7'); 
O| 为 实测 均值 (ms ) ;n 为 实测 数据 量 "。 
2.3 径流 突变 检验 

变化 环境 下 ,乌拉 盖 河 流域 水 文 要 素 无 疑 会 受 
到 一 定 程度 的 干预 和 扰动 ,本 研究 采用 Mann-Kend- 
al 突变 检验 法 用 于 径流 突变 分 析 。 与 此 同时 ,应 用 
降水 -径流 量 双 累 积 曲线 作 为 突变 点 的 辅助 检验 。 
Mann-Kendall 趋势 检验 法 (简称 M-K 法 ) 是 由 Mann 
和 Kendall 提出 并 改进 的 统计 方法 ,经 常用 于 分 析 降 
水 径流 和 气温 等 要 素 的 长 时 间 序 列 变化 "*"。 其 特 
点 在 于 对 样本 数据 的 分 布 无 要 求 , 受 异常 值 的 干扰 
较 小 ,在 水 文 . 气 象 等 数据 中 适用 性 良好 , 且 计 算 简 
便 。 当 M-K 检 验 突变 点 时 ,检验 变量 5$ 计 算 公 式 
如 下 : 


R= 


gess i (k=2,3,4, +++, n) (4) 
i J 


azio ari. gejai (5) 
loner 
UB,=-UF, ， 

23, 6 
人 ( i W 


式 中 :秩序 S, 是 第 i 时 刻 数值 大 于 j 时 刻 数值 个 数 
的 累计 值 。 通 过 分 析 曲 线 UB, 和 UF 可 以 表明 径 
流 深 的 变化 趋势 又 可 以 明确 其 发 生 的 突变 年 份 。 
UF , 值 的 正 负 代表 径流 深 的 增加 或 减少 趋势 , UB 
和 UF 是 否 越过 临界 值 线 则 表示 变化 趋势 的 显著 
GA, UB, 和 UF 的 交点 位 于 临界 线 内 ,该 交点 
为 突变 点 号 。 

降水 -径流 量 双 累积 曲线 是 检验 两 个 参数 间 关 
系 一 致 性 及 其 变化 的 常用 方法 。 所 谓 双 累积 曲线 
就 是 在 直角 坐标 系 中 绘制 的 同期 内 一 个 变量 的 连 
续 味 积 值 与 另 一 个 变量 连续 累积 值 的 关系 线 , 可 用 
于 水 文 气象 要 素 一 致 性 的 检验 、 缺 值 的 插 补 或 资料 
校正 ,以 及 水 文 气象 要 素 的 趋势 性 变化 及 其 强度 的 
分 析 '"。 本 研究 以 降水 和 径流 深 为 变量 作 双 累积 
曲线 ,将 曲线 斜率 发 生 明显 变化 的 年 份 定义 为 突变 
年 ,突变 年 以 前 作为 基准 期 ,也 称 人 类 活动 稳定 期 ; 
突变 年 以 后 作为 突变 期 , 称 人 类 活动 剧烈 期 。 
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2.4 区 分 气候 变化 与 人 类 活动 对 径流 变化 的 影响 
为 了 解析 近 40 a 气候 变化 与 人 类 活动 对 乌拉 六 
河流 域 的 综合 影响 ,如 暧 干 化 趋势 及 过 度 放牧 造成 
的 草地 退化 。 本 研究 选取 乌拉 盖 河 流域 1981 一 
2020 年 的 气象 数据 以 及 过 去 5 个 时 期 (1980 年 、 
1990 年 .2000 年 .2010 年 和 2020 年 ) 土 地 利用 / 禾 被 ， 
结合 情景 设置 法 模拟 不 同 变化 环境 下 的 径流 量 。 
共 设 置 11 种 情景 ,分 别 用 sl ~ s11 表 示 ( 表 1)。 人 情景 
1(4s1) 将 基准 期 的 气象 数据 和 土地 利用 / 履 被 驱动 模 
型 ,得 到 基准 期 土地 利用 / 覆 被 和 气候 条 件 下 的 径流 
模拟 值 。 为 保证 突变 期 模拟 值 和 基准 期 的 径流 量 
在 物理 成 因 上 的 一 致 性 ,情景 2 和 情景 3 将 保持 率 
定好 的 模型 参数 不 变 , 模 型 分 别 输入 突变 期 气象 数 
据 和 土地 利用 / 窗 被 ,模拟 得 到 突变 期 气候 和 土地 利 
用 / 覆 被 条 件 下 的 径流 值 。 为 了 进一步 了 解 不 同时 


总 变化 值 ; W, 气候 变化 对 径流 量 的 贡献 率 ; W, 人 类 
活动 对 径流 量 的 贡献 率 。 


3 结果 与 分 析 


3.1 SWAT 模 型 校准 与 验证 

本 研究 建立 的 SWAT 模 型 将 乌拉 盖 河 流域 划分 
成 41 个 子 流 域 和 153 个 水 文 响应 单元 。SWAT 模 型 
参数 较 多 ,采用 SWAT-CUP 软件 对 参数 进行 敏感 性 
分 析 , 选 取 最 敏感 的 14 个 参数 进行 重点 率 定 ,得 到 
的 敏感 参数 校准 值 见 表 2。 模 型 预 热 期 为 1976 一 
1980 年 , 率 定 期 为 1981 一 1998 年 ,验证 期 为 1999 一 
2012 年 。 由 站 点 实测 径流 数据 与 模拟 径流 数据 对 
比 来 看 ,径流 模拟 率 定 期 和 验证 期 评价 指标 NSE 和 
尼 均 高 于 0.6,PBLAS 小 于 20%( 图 2a)。 乌 拉 盖 河 是 


段 气 候 变化 和 人 类 活动 对 径流 影响 的 差异 ,本 研究 
将 乌拉 盖 河 流域 过 去 40 a 的 气候 和 十 地 利用 /有 覆 被 
分 为 4 个 时 段 (1980 一 1990 年 、1990 一 2000 年 、 
2000—2010 年 .2010 一 2020 年 ) ,将 其 设置 为 s4~ 
s11 等 8 个 情景 。 通 过 以 上 11 种 情景 分 别 计算 气候 
以 及 土地 利用 / 黎 被 引起 的 径流 变化 值 ,具体 表达 
如 下 : 


AQ=AQ,+AQ, (7) 
_ AQ. 
oar (8) 
_ AQ, 
= Ro (9) 


AP: AQ, 为 基准 期 气候 变化 引起 的 径流 变化 值 ; 
0, 为 突变 期 人 类 活动 引起 的 径流 变化 值 ; AQ 径流 


季节 性 河流 ,实测 径流 表现 为 夏季 径流 量 较 大 , 冬 
季 断 流 的 特点 。 改 进 的 SWAT 模 型 能 够 较 好 的 捕捉 
径流 的 季节 性 分 布 特征 。 此 外 ,通过 对 比 MODIS 数 
据 与 模型 模拟 的 积 雪 面 积 . 蒸 散发 结果 发 现 ,模拟 
值 与 遥感 数据 拟 合 较 好 , 且 两 者 的 相关 系数 均 高 于 
0.8( 图 2b、 图 2c)。 根 据 多 目标 评定 结果 来 看 ,模型 
模拟 结果 较 好 ,可 用 于 下 一 步 研究 。 
3.2 乌拉 盖 河 流域 径流 变化 时 空 特征 

根据 SWAT 模 拟 结果 ,1981 一 2020 年 乌拉 盖 河 
流域 年 均 径 流 深 以 1.19 mm- a` HIE K E aD 
趋势 (图 3a)。 其 中 ,径流 深 的 最 大 值 出 现在 1998 年 
(183.24 mm) , 最 小 值 出 现在 2017 年 (12.08 mm) 。 
不 同 河 段 径流 变化 趋势 与 整个 流域 变化 一 致 ,径流 
最 大 值 均 出 现在 1998 年 ,但 径流 大 小 相差 甚大 (图 


表 1 模型 模拟 情景 设置 


Tab.1 Model simulation scenario settings 


情景 气象 数据 土地 利用 / 巷 被 数据 目的 

sl 1980 一 2000 年 1980 年 计算 乌拉 盖 河 流域 突变 年 前 后 基准 期 和 突变 期 气候 变化 与 人 类 活动 对 径流 贡献 率 
s2 2000—2020 年 1980 年 

s3 2000 一 2020 年 2020 年 

s4 1980 一 1990 年 1980 年 计算 1980 一 1990 年 气候 变化 与 人 类 活动 对 径流 的 贡献 率 
s5 1980—1990 íE 1990 4F 

s6 1990—2000 4 1990 4E 计算 1990 一 2000 年 气候 变化 与 人 类 活动 对 径流 的 贡献 率 
s7 1990 一 2000 年 2000 年 

s8 2000 一 2010 年 2000 年 计算 2000 一 2010 年 气候 变化 与 人 类 活动 对 径流 的 贡献 率 
s9 2000 一 2010 年 2010 年 

s10 2010—2020 年 2010 4E 计算 2010 一 2020 年 气候 变化 与 人 类 活动 对 径流 的 贡献 率 
sll 2010—2020 4 2020 年 
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表 2 参数 敏感 性 分 析 结 果 


Tab.2 Parameter sensitivity analysis results 


序号 参数 名 称 参数 含义 取 值 范围 最 终 取 值 
1 CN2 SCS 径 流 曲 线 系数 30~98 85.19 
2 ALPHA_BF 基 流 消退 系数 /d 0~1 0.11 
3 GW_DELAY 地 下 水 延迟 系数 /d 0~500 4.36 
4 GWQMN 浅 地 下 水 产 流 阔 值 /mm 0~1000 104.25 
5 SFTMP 降雪 温度 /%C -10~10 -1.03 
6 SMTMP 融雪 基准 温度 /%C -10~10 2.83 
7 SMFMX 最 大 积 雪 融化 速率 /(mm%C7') -20~20 8.29 
8 SMFMN 最 小 积 雪 融化 速率 /(mm" C) -20~20 2.26 
9 SNOCOVMX 100% 积 雪 禾 盖 时 所 对 应 的 积 雪 含水 量 /mm 0~500 161.40 
10 ESCO 土壤 层 蒸 发 补偿 系数 0~1 0.65 
11 EPCO 植物 吸收 补偿 因子 0~1 0.11 
12 SOL_K 饱和 水 力 传导 系数 /mm'h ) 0~2000 57.90 
13 SOL_AWC 土壤 有 效 持 水 量 /mm 0~1 0.36 
14 SURLAG Hee iin at AY Ti] 0~24 13.79 
(a) 径流 量 模拟 对 比 
< > 率 定 期 < >< > 验证 期 < > 
140. i 
一 一 模拟 数据 率 定期 : 验证 期 : 
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图 2 乌拉 羡 河 流域 SWAT 模 拟 效果 
Fig. 2 SWAT simulation effect of the Ulagai River Basin 


3b)。 通 过 M-K 检 验 结果 发 现 ,UF 和 UB 两 条 曲线 
在 2000 年 出 现 交点 ,日 通过 显著 性 检验 (图 3c), 表 
明 在 2000 年 乌拉 盖 河 流域 径流 量 发 生 明 显 突变 。 
基于 降水 -径流 量 双 累积 曲线 法 辅助 验证 突变 年 ， 


在 2000 年 以 后 曲线 的 曲率 出 现 明显 偏 移 (图 3d), 同 
样 表明 径流 深 在 2000 年 发 生 突变 。 这 与 研究 区 附 
近 其 他 河流 的 研究 结果 一 致 “”。 

通过 分 析 乌 拉 盖 河流 域 多 年 平均 径流 深 的 空 
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图 3 乌拉 羡 河 流域 1981 一 2020 年 径流 时 间 变 化 以 及 突变 检验 
Fig.3 Runoff time change and mutation test in the Ulagai River Basin from 1981 to 2020 


间 分 布 情况 ,发 现 流域 中 部 径流 较 大 ,东部 径流 较 
小 ,并 且 由 东 向 西 呈 先 增 大 后 减 小 的 分 布 格局 , 流 
域 上 .中游 径流 较 小 ,下 游 径流 较 大 (图 4a)。 由 马 
拉 盖 河流 域 径流 深 变 化 的 空间 分 布 图 来 看 (图 4b)， 
整个 流域 径流 深 均 呈 不 同 程度 的 增加 趋势 , 且 变 化 
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0~15 
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fim 50~120 
ME 120~520 
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幅度 大 小 与 时 空 分布 格 局 相近 ,下 游 低 海拔 区 域 较 
上 中 游 增加 趋势 更 为 明显 。 

3.3 乌拉 盖 河 流域 径流 归 因 分 析 

3.3.1 乌拉 盖 流 域 不 同 河 段 径流 变化 归 因 分 析 根 
据 突 变 检 验 结 果 将 研究 时 段 分 为 基准 期 (1981 一 
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4 乌拉 羡 河 流域 1981 一 2020 年 径流 深 及 变化 的 空间 分 布 
Fig.4 Spatial distribution of runoff depth in the Ulagai River Basin from 1981 to 2020 and spatial distribution of variation 
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1999 年 ) 和 突变 期 (2000 一 2020 年 ) 。 再 结合 情景 
sl~s3 结果 发 现 , 气 候 变 化 使 年 均 径 流量 减少 40.07 
mm ,贡献 率 为 93.84%; 人 类 活动 使 年 均 径 流量 减少 
1.74 mm, 贡献 率 4.16%。 因 此 ,气候 变化 是 乌拉 盖 
河流 域 径 流 变 化 的 主要 原因 。 乌 拉 盖 河流 域 自 基 
准 期 之 后 气温 持续 增加 、 降 水 量 减少 ,径流 量 也 随 
之 减少 ,突变 期 较 基 准 期 的 降水 量 减 少 了 14% , 气 
温 增 加 了 0.82 % ,总 径流 量 减少 了 69.20%。 降 水 量 
减少 .气温 增加 使 径流 量 减 少 ,但 是 多 年 径流 深 的 
相对 变化 远大 于 降水 量 和 气温 的 变化 。 这 表明 人 
类 活动 对 径流 量 的 影响 程度 也 不 容 忽 视 ”。 如 图 5 
所 示 ,气候 变化 和 人 类 活动 对 于 乌拉 盖 河 流域 不 同 
河 段 的 影响 不 尽 相 同 , 人 类 活动 对 上 中、 下 游 的 贡 
献 率 分 别 为 1.69% 、4.36% 和 5.03%。 


m 气候 变化 o 人 类 活动 


1007 95.84 38.31 95.64 94.97 


80} 


20} 


整个 流域 上 游 中 游 下 游 


图 5 乌拉 盖 河 流域 径流 变化 归 因 
Fig. 5 Attribution of runoff changes in the Ulagai River Basin 


3.3.2 乌拉 盖 河 流域 不 同时 段 径 流 变化 归 因 分 析 为 
了 明晰 不 同时 间 段 内 气候 变化 和 人 类 活动 对 乌拉 善 
河 径流 的 影响 程度 及 差异 ,根据 5 期 土地 利用 / 履 被 
数据 将 研究 时 期 划分 为 不 同 的 4 个 时 段 。 结 合 情 景 
(s4~s11) 分 别 研 究 1980—1990 年 1990—2000 年 、 
2000 一 2010 年 和 2010 一 2020 年 不 同时 段 内 的 径流 
演变 特征 ,量化 气候 变化 和 人 类 活动 对 径流 变化 的 
贡献 。 可 以 看 出 ,不同 时段 气候 变化 和 人 类 活动 对 
流域 径流 的 贡献 率 均 不 相同 ,并且 有 正 负 两 种 影响 
(图 6)。 在 4 个 不 同时 段 内 气候 变化 对 径流 量 的 贡 
献 率 逐 渐 减 少 ,由 88.26% 减 少 到 25.47%。 而 在 
2000 一 2010 年 内 气候 变化 和 人 类 活动 均 对 径流 产 
生 了 负面 影响 ,贡献 率 分 别 为 -48.41% 和 -51.19%。 
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图 6 乌拉 盖 河 流域 不 同时 段 气候 变化 与 人 类 活动 对 径流 
的 贡献 率 


Fig. 6 Contribution of climate change and human activities 


to runoff at different times in the Ulagai River Basin 


4 讨论 

对 于 水 资源 缺乏 的 干旱 半 干 旱 牧 区 草原 内 陆 
河流 域 ,降水 成 为 了 主要 水 量 补给 来 源 。 总 体 而 
言 ,1981 一 2020 年 乌拉 盖 河 流域 降水 量 平 均值 为 
323.89 mm, 在 1998 年 达到 峰值 (552.73 mm) , 以 
8.91 mm* (10a) "的 速率 旺 减 少 趋 热 。 近 40 a 气温 以 
波动 形式 上 升 ,升温 速率 可 达 0.43 % .(10a) ,使 该 
流域 呈 暧 干 化 趋势 (图 7)。 研 究 显示 ,流域 暧 干 化 
的 背景 下 径流 量 以 11.93 mm (10a)“! 的 速率 减少 。 
而 M-K 检 验 和 降水 -径流 量 双 累积 曲线 的 结果 均 显 
示 年 均 径 流 深 在 2000 年 发 生 突变 ,这 表明 乌拉 羡 河 
流域 径流 量 在 2000 年 之 前 主要 受气 候 变 化 的 影响 ， 
而 后 受气 候 和 人 类 活动 的 共同 影响 。 流 域 附近 的 
其 他 流域 径流 也 存在 着 突变 现象 ,日 突变 期 总 径流 
量 小 于 基准 期 总 径流 量 "*”。 表 示 突 变 年 后 径流 量 
的 减少 受气 候 变 化 及 过 度 放牧 .盲目 开垦 以 及 水 利 
水 库 建设 等 人 类 活动 的 双重 影响 。 乌 拉 羡 河流 域 
暧 干 化 的 趋势 ,在 一 定 程度 上 导致 径流 量 减 少 。 然 
而 ,突变 期 较 基 准 期 多 年 径流 深 的 相对 变化 远大 于 
降水 量 和 气温 的 变化 ,凸显 了 过 度 放牧 .盲目 开采 、 
以 及 不 合理 的 农业 灌溉 等 人 类 活动 对 径流 的 负面 
影响 。 在 经 济 快速 发 展 和 人 口 持续 增加 的 过 程 中 
城市 扩张 .工矿 用 地 增多 .道路 建设 均 会 产生 耗 水 
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图 7 乌拉 盖 河 流域 气候 要 素 及 径流 变化 


Fig. 7 Changes in climatic elements and runoff runoff in the Ulagai River Basin 


SIM! 。 据 统计 ,乌拉 盖 河 流域 建设 用 地 面积 
1980 年 的 53.66 km 增加 到 2020 年 的 264.13 km’. 
此 外 ,该 流域 畜牧 业 更 加 发 达 , 人 研究 时 段 内 牲畜 数 
量 持 续 增加 ,2000 年 以 4x10' 头 的 数量 达到 近 40 a 
最 高 值 。 过 度 放牧 等 不 合理 的 人 类 活动 导致 牧区 
草原 植被 退化 ,使 雨水 拦蓄 能 力 下 降 , 从 而 加 速水 
土 流失 i。 而 人 类 活动 对 上 中、 下 游 的 影响 也 不 
尽 相 同 , 贡 献 率 分 别 为 1.69% 、4.36% 和 5.03%。 不 
同 河 段 的 土地 利用 / 覆 被 类 型 工业 以 及 畜牧 业 发 展 
的 不 同 程度 可 能 会 导致 这 种 差异 。 如 流域 上 游 长 
期 以 来 以 农业 发 展 为 主 , 且 降 水 量 较 高 温度 较 低 ， 
径流 受 人 类 活动 的 影响 较 小 。 而 以 草地 覆 被 为 主 
的 下 游 区 域 , 主 要 依靠 畜牧 业 为 生 。 而 下 游 性 畜 数 
量 的 增加 以 及 建立 工矿 等 人 类 生活 方式 使 耗 水 量 
增加 ,从 而 减少 水 资源 1。 

通过 模拟 乌拉 盖 河 流域 5 个 时 期 不 同 土地 利 
用 / 覆 被 条 件 下 的 径流 过 程 并 结合 下 垫 面 变化 以 及 
政府 实施 的 工程 措施 ,5 个 时 段 的 产 流 比如 表 3 所 
示 ,1980 一 2020 年 研究 区 暧 干 化 趋势 下 流域 径流 呈 
先 增加 后 减少 趋势 ,1981 一 2000 年 径流 均 呈 增加 趋 


势 ,增加 了 18.45 mm;2000—2020 年 呈 减 少 趋势 ,总 
径流 减少 了 39.85 mm。 上 述 分 析 表 明 ,突变 年 
(2000 年 ) 后 人 类 活动 是 乌拉 盖 河 流域 径流 显著 减 
少 的 重要 因素 ,而 土地 利用 / 履 被 变化 是 人 类 活动 的 
集中 体现 形式 之 一 *。 

由 乌拉 羡 河 流域 近 40 a 土地 利用 / 覆 被 可 知 (图 
8) ,发 现 1981 一 2020 年 流域 内 耕地 面积 先 增 后 减 ， 
林地 和 草地 面积 变化 趋势 与 之 相反 ,这 与 乌拉 盖 河 
流域 毁 林 毁 草 开 况 到 退耕 还 林 还 草 相 对 应 |。 
1981 一 2000 年 径流 持续 的 增加 ,其 原因 为 该 时 段 流 
域 发 展 第 一 产业 ,流域 牲畜 总 量 为 2000 一 2020 年 的 
PER., TEF E 超载 养 冀 以 及 灌溉 量 的 增 
加 等 导致 下 热 面 的 改变 ,一 定 程度 上 削弱 了 其 对 雨 
水 的 拦蓄 作用 ,使 径流 量 增加 '"。 此 外 ,突变 分 析 
结果 显示 ,与 基准 期 相 比 降幅 达到 69.20%。 尽 管 
2000 年 开始 政府 采取 一 系列 水 土 保持 措施 ,导致 耕 
地 面积 减少 以 及 玻 林地 和 低 履 盖 度 草地 面积 的 增 
加 ,在 一 定 程度 上 减少 了 出 口 断面 的 径流 量 ,防止 
了 水 土 流 失 中 。 然 而 ,同时 间 内 锡林郭勒 盟 第 二 、 
三 产业 迅猛 发 展 ,建立 乌拉 盖 水 库 为 附近 的 煤矿 、 


表 3 乌拉 盖 河 流域 不 同时 段 产 流 比 


Tab.3 Ratio ofproduction to flow in the Ulagai River Basin at different times 


年 份 人 类 活动 方式 地 表 径 流 /mm 二 中 流 /mm 地 下 径流 /mm 总 产 水 量 /mm 
1980 年 原始 状态 8.41 18.12 18.02 45.55 
1990 年 过 度 开垦 20.69 15.28 14.26 50.91 
2000 年 过 度 放牧 12.60 31.82 23.00 64.00 
2010 年 退耕 还 林 还 草 3.28 11.80 1.10 15.19 
20204 持续 水 土 保持 工程 8.30 13.34 2.40 24.15 
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图 8 乌拉 盖 河 流域 土地 利用 / 窗 被 变化 


Fig. 8 Land use/cover change in the Ulagai River Basin 


石油 的 开采 提供 工业 用 水 2 ,年 度 工业 用 水 量 达到 
T 4750 m。 这 极 大 地 消耗 了 乌拉 羡 河 流量 ,导致 
2010 年 后 径流 量 急速 下 降 。 因 此 ,径流 的 减少 不 仅 
受 土地 利用 /和 覆 被 ,还 包括 人 口 增 长 和 经 济 发 展 导致 
对 水 资源 的 开发 利用 等 因素 。 流 域内 从 2000 “EFF 

台 实施 了 一 系列 生态 保护 措施 ,包括 “退耕 还 林 还 
草 ”“ 京 津 风 沙 源 治理 ”,“ 未 封禁 牧 ” 等 水 土 保持 措 
We?” ,但 径流 量 在 2000 年 依旧 呈 增 加 趋势 , 而 在 
2010 年 开始 减少 ,这 表明 水 土 保持 措施 具有 一 定 的 
滞后 效应 。 


5 结论 


(1) SWAT 模 型 在 乌拉 盖 河 流域 的 径流 模拟 适 
用 性 良好 , 率 定期 及 验证 期 径流 模拟 结果 与 实测 径 
流 的 NSE 及 尼 均 在 0.62 以 上 ,PBLAS 小 于 18.8% 。 
除 此 之 外 ,MODIS 积 雪 面积 . 蒸 散 发 与 模型 模拟 值 
的 拟 合 系数 均 高 于 0.80。 

(2) XE 40 a 降水 减少 及 气温 增加 的 背景 下 , 乌 
拉 立 河流 域 年 均 径 流 深 以 1.19 mm: a! AYR SE 


著 减 少 趋 势 , 晶 在 2000 年 发 生 突 变 。 空 间 上 由 东 向 
西 呈 先 增 大 后 减 小 的 分 布 格局 。 

(3) 乌拉 盖 河 流域 气候 变化 是 径流 变化 的 主要 
驱动 因子 ,对 径流 的 贡献 率 为 95.84%。 对 不 同 河 段 
的 贡献 率 也 有 所 差异 ,对 上 中、 下 游 的 贡献 率 分 别 
为 98.31% .95.64% Fil 94.97% 。 对 于 不 同时 间 段 , 流 
域 气候 变化 对 径流 的 贡献 率 逐 渐 减 少 , 由 1980 年 的 
88.26% 减 少 到 2020 年 的 25.47%。 而 不 同人 类 行为 
活动 导致 径流 变化 趋势 及 幅度 也 有 较 大 的 差异 。 
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Abstract: The Ulagai River Basin in Inner Mongolia is a typical inland river basin in the pastoralist grasslands of 
the arid and semi-arid regions with extremely fragile ecosystems, and climate warming and human activities can 
directly affect hydrological changes. In this study, we used the improved SWAT (Soil and Water Assessment 
Tool) hydrological model, M-K trend test, precipitation-runoff double accumulation curve, and scenario analysis 
to systematically analyze the spatial and temporal characteristics of runoff in the Ulagai River Basin from 1981 to 
2020 and to quantify the differences between climate change and human activities on runoff in different periods 
and river sections. The differences in the impact of climate change and human activities on the runoff in different 
periods and river sections were quantified. The results show that the SWAT model significantly impacted runoff 
in the Ulagai River Basin from 1981 to 2020. The results show that the SWAT model has good applicability in the 
Ulagai River Basin, with NSE and R’ above 0.62 and PBLAS less than 18.8% for both the periodic and validation 
periods. Under the warm and dry trend of the basin in the past 40 years, the runoff volume decreased significantly 
in the upper, middle, and lower reaches, then a sudden change occurred in 2000. Considering the contribution of 
human activities such as climate change and overgrazing, blind reclamation and construction of water 
conservancy reservoirs to the change of runoff in the basin was 95.84% and 4.16%, respectively. Also, the 
contribution of human activities to different river segments in the basin varied, with 1.69%, 4.36%, and 5.03% 
from upstream to downstream. The contribution rates also differed significantly across different periods, from 
88.26% in 1980 to 25.47% in 2020. The trend and magnitude of runoff changes due to different human activity 
patterns in different periods also differed. This study’ s results can provide a reference basis for the sustainable 
use and rational scheduling of water resources in the inland river basin of pastoral grasslands. 


Keywords: runoff; climate change; human activities; SWAT model; Ulagai River Basin 


